The Eucalyptus grandis chloroplast proteome: Seasonal variations in leaf development

Chloroplast metabolism is very sensitive to environmental fluctuations and is intimately related to plant leaf development. Characterization of the chloroplast proteome dynamics can contribute to a better understanding on plant adaptation to different climate scenarios and leaf development processes. Herein, we carried out a discovery-driven analysis of the Eucalyptus grandis chloroplast proteome during leaf maturation and throughout different seasons of the year. The chloroplast proteome from young leaves differed the most from all assessed samples. Most upregulated proteins identified in mature and young leaves were those related to catabolic-redox signaling and biogenesis processes, respectively. Seasonal dynamics revealed unique proteome features in the fall and spring periods. The most abundant chloroplast protein in humid (wet) seasons (spring and summer) was a small subunit of RuBisCO, while in the dry periods (fall and winter) the proteins that showed the most pronounced accumulation were associated with photo-oxidative damage, Calvin cycle, shikimate pathway, and detoxification. Our investigation of the chloroplast proteome dynamics during leaf development revealed significant alterations in relation to the maturation event. Our findings also suggest that transition seasons induced the most pronounced chloroplast proteome changes over the year. This study contributes to a more comprehensive understanding on the subcellular mechanisms that lead to plant leaf adaptation and ultimately gives more insights into Eucalyptus grandis phenology.


Introduction
Chloroplast is a subcellular compartment in which solar radiation is used to generate sugar compounds in a complex process named photosynthesis. It is also hosts several other metabolic processes that are essential for plant survival and is composed of inner and outer membrane, stroma, thylakoidal membrane, and lumen [1]. Moreover, Eucalyptus leaves have a great amount of antioxidants, particularly flavonoids, which provide protection against from ROS damage, oxidative stress and can increase stress tolerance [2]. Plastid genome encodes 120-130 proteins, however approximately 2000 proteins, are synthesized in the cytosol, and imported into the organelle through the envelope membranes [3,4]. In Eucalyptus grandis, the chloroplast genome presents approximately 160 kb [5,6].

Fluorescence measurement and chlorophyll relative quantification
Fluorescence analysis was performed using the fluorometer FluorPen FP 100 (Photon Systems Instruments) assessing the adaxial surface of leaves from different regions of the 15 branches.
The fluorescence values obtained directly from the device were used to assess the potential quantum efficiency of photosystem II (FV/FM). Relative chlorophyll was quantified on CCM-200 equipment (Opti-Sciences), and the adaxial surface of leaves from different regions of the 15 branches was examined. The CCI values obtained directly from the device were used to assess relative chlorophyll quantification. Statistical analysis was performed according to the ANOVA test (p < 0.05), and Tukey's test was used for comparison among different leaves.

Chloroplast proteome extraction
Six grams of leaves were collected, fine sliced, and transferred to tubes containing 20 mL of chloroplast isolation buffer (CIB, CP-ISO-Sigma Aldrich). The leaf tissue was homogenized (Omni GLH Turrax) and filtered through a nylon membrane. The processed material was centrifuged, and the chloroplasts were isolated in an 80%/40% (1:2 v/v) Percoll density gradient. After washing, the chloroplast pellet was resuspended in a protein extraction buffer (Thiourea 2M; CHAPS; Tris-HCl 1M, pH 8; Glycerol 50%; Triton X-100; DTT 1M and Urea 8M) and the lysis of intact chloroplasts was carried out through 20 cycles of sonication. After centrifugation, the supernatant was precipitated in cold acetone for 16 hours at -30˚C, and the protein pellet obtained was solubilized in Tris buffer (125 mM Tris pH 6.8, 20% glycerol, 1% SDS, 0.1% DTT).

Total leaf proteome extraction
To assess the enrichment of the chloroplast subproteome in relation to the total leaf proteome, intact leaves were frozen in -80˚C and grinded in Tris extraction buffer (125 mM Tris, pH 6.8, 20% glycerol, 1% SDS, 0.1% DTT). Upon centrifugation at 4˚C for 15 min, the supernatant was transferred to a clean tube and the proteins were precipitated in cold acetone for 16 hours at -30˚C. Finally, the proteins were solubilized in Tris buffer before protein quantification assay.

Protein quantification assay and SDS-PAGE
Sample quantifications followed the Bradford method using a BSA curve as reference. Aliquots of 50 μg of proteins were separated for 1 hour by denaturing electrophoresis on a polyacrylamide gel containing sodium dodecyl sulfate (SDS-PAGE) to obtain a single band. The electrophoretic system was composed of a stacking gel (5%) and a separation gel (12.5%). The single protein band was excised from the gel and subjected to in-gel digestion.

Mass spectrometry analysis
Proteins were in-gel digested using trypsin at 1:10 proportion (trypsin:protein) [36]. To clean up the peptide samples, extracts were filtered in an acetate cellulose membrane (Spin-X, Corning) by centrifugation for 7 minutes. The flow-through was collected using a clean tube, concentrated, and then desalted with ZipTip Pipette Tips (Millipore) according to the manufacturer's protocol. The amount of 1ug of extracted peptides was analyzed through liquid chromatography coupled to a tandem mass spectrometer. Tryptic peptides were resuspended in H 2 O containing 0.1% formic acid and separated in an EASY-nLC1000 chromatographic system (Thermo Fisher Scientific) for 95 minutes using a C18 nano-column under a constant flow of 400 nL/min and increasing concentrations of acetonitrile (7 to 30% over 70 min and then 70% in 25 min). The peptides were ionized through electrospray (3.5 kV) and the mass spectrometer (Q-Exactive, Thermo Fisher Scientific) was operated in data-dependent acquisition mode according to a top 20 acquisition method with dynamic exclusion time adjusted to five seconds. Full scans and MS/MS scans were acquired with resolutions equal to 70,000 and 35,000 FWHM (full width at half maximum), respectively. The MS/MS spectra were generated from HCD fragmentation (high-energy collision dissociation) of the precursor peptide ions isolated under 35 eV collision energy. Mass spectrometric data are available via ProteomeXchange under identifier PXD029004.

Protein identification
The peptides and protein inferences were identified through stringent searches of experimental data against the protein database Eucalyptus grandis v.2.0 (http://www.phytozome.net/) on MaxQuant (v. 1.6.3.3) [37] based on dedicated search parameters: 20 ppm tolerance for precursor ions, oxidation of methionines and carbamidomethylation of cysteines-selected as dynamic and static modifications, respectively. Maximum missed cleavages of 2 and false-discovery rate adjusted to 1%. Relative abundance of the chloroplast proteins were quantified based on the label free-based approach (LFQ intensity). Relative abundance of total leaf proteins was obtained following the NSAF ("normalized spectral abundance factor") approach [38]. Mass spectrometry proteomics data are currently stored on the ProteomeXchange Consortium via PRIDE partner repository under dataset identifier PXD029004.

Statistical analysis and functional annotation
Data processing and statistical analysis were carried out on the Perseus software (v. 1.6.15.0) [39]. Firstly, all contaminants and reverse sequences were removed from data analysis. Subsequently, label-free quantification data (LFQ) were transformed into log 2 scale, and a protein was considered a true hit if identified in at least two out of three replicates. Missing values imputation was carried out from the normal distribution curve from each biological replicate. The means of ANOVA test and Tukey comparison test (p <0.05) were used in the comparative analyses. Multivariate analyses were carried out with differentially abundant proteins for both assays on the Perseus software. For PCA and hierarchical clustering, data were z-score normalized. Hierarchical clustering was performed using Euclidean distance and complete linkage. Functional annotation of the differentially abundant proteins was carried out using the Blast2GO platform (p< 0,05) [40]. The ten most representative level-four biological processes were considered in the functional characterization aiming to balance between a comprehensive overview and detailed information on the targeted biological process. Level four was chosen for including photosynthesis, an important process in the context of leaves and chloroplasts.

Subcellular prediction
In silico analysis was performed using predictors of subcellular location and chloroplast proteome repositories to assess and confirm the cellular location of the identified proteins. For this purpose, a protein was considered plastidial if identified by at least one of the subcellular prediction algorithms-ChloroP [41] and Predotar [42]-or by presenting at least one homologous protein in any of the databases used-PPDB [43] or AT-Chloro [44].

Results and discussion
Before performing any description or comparative proteome analysis, we decided to verify the efficiency of the chloroplast enrichment strategy used herein. For this purpose, we compared the abundance levels, in terms of NSAF and LFQ intensities, of a well-known chloroplast proteoform Eucgr.C03525.1 (Ribulose-bisphosphate carboxylase-RBCL) against five potential mitochondrial contaminants commonly found in chloroplast extracts (Table 1). Quantitative analysis indicated that the chloroplast proteoform was enriched upon the extraction process, whereas the mitochondrial proteoforms presented a decrease in the abundance levels. These results indicated that our chloroplast extraction and enrichment strategy was effective in reducing the cross-contaminants protein abundance, as well as in the enrichment process of chloroplast-related proteins in the target samples.

Photosynthetic proteins comprise most of the Eucalyptus grandis chloroplast proteome
Aiming to illustrate a general overview of the chloroplast proteome isolated from Eucalyptus grandis leaves, we compiled the identification datasets from both leaf development and seasonal experiments. The combined analysis resulted in a stringent, non-redundant dataset of 431 chloroplast proteins (S2 Table). Gene ontology analysis indicated that most of the isolated proteins were involved in the "organonitrogen compound metabolic process", "organic substance biosynthetic process", and "macromolecule metabolic process" (Fig 1A). Although these biological processes were the most prominent in terms of protein identifications, an MS1-based quantitative analysis indicated that the "photosynthesis process" represented 11% of the isolated chloroplast proteomes (Fig 1B). The fraction of "photosynthetic process" related to proteins becomes even more pronounced within the isolated proteome upon the analysis of the most abundant proteins identified in the entire experiment ( Fig 1C). The RuBisCO small chain (Eucgr.J01502.2.p) was the most abundant protein identified herein, showing approximately 4000 MS/MS counts. The proteins involved in the photosystem II structure were also highly abundant: two extrinsic proteins (PsbO-2 -Eucgr.I01025.1.p and PsbQ-2 -Eucgr. D00854.1.p) that form the oxygen-evolving complex (OEC) and a light-harvesting complex protein from the PSII (Eucgr.F04099.1.p) corresponding to 1500 to 3000 MS/MS counts ( Fig 1C).
Regarding the suborganellar localization of the 431 chloroplast proteins, 137 proteins were present in the stroma compartment, whilst 135 were thylakoidal proteins. A total of 87 proteins were predicted in silico to be located in the chloroplast outer membrane, and 13 identifications were found to be part of the chloroplast inner membrane. Finally, 52 proteins were located in the chloroplast lumen and 7 proteins herein identified had been previously reported as part of the plastoglobule structure ( Fig 1D).

Chloroplast proteome reflects the physiological status in developing leaves
Based on the varying FV/FM ratios and the differences observed from the relative chlorophyll quantitation analysis (S2 Fig), we isolated the chloroplast proteome from young, middle, and mature leaves. Quantitative data from the isolated proteins were used to demonstrate the changes within this subcellular structure throughout the leaf development process. Pairwise analysis indicated that 61 proteins were considered to be differentially abundant over the development process, with the young chloroplast proteome as the most dissimilar of the three different regions (Fig 2A). Such result is in accordance with physiological evidence highlighted in S1 Fig, as the young leaf tissue presented the most significant variations concerning middle and mature regions. A study on sugarcane leaves have also demonstrated that physiological parameters, such as chlorophyll content, nitrogen content, and RubisCO activity, were greater and more consistent between mature and transition tissues than in young regions [45]. Clustering analysis allowed to classify the differentially abundant proteins (DAP) in three distinctive groups: upregulated-17 proteins with higher abundance over leaf development, downregulated-a set of 29 proteins with lower abundance along the studied process, and an undefined group, comprising 15 proteins with an unclear abundance pattern ( Fig 2B).
Aiming to highlight changes in the chloroplast metabolic processes along leaf development, up and downregulated proteins were cross-referenced to corresponding gene ontology terms. Functional categories revealed differences between the chloroplast protein groups (Fig 2C). Leaf development induced "cellular catabolic processes" in Eucalyptus grandis chloroplasts for being the most representative category within the upregulated group. Responses to stress, mainly to oxidative stress and cellular detoxification, were the most abundant functional category in upregulated throughout leaf development, based on the maturity of the photosynthetic apparatus and the higher photosynthesis rates in young leaves. In non-stressing conditions, mature leaves are expected to maintain cellular homeostasis with a dynamic balance between the light energy absorbed and energy consumption [46]. In mature leaves, generation of reactive oxygen species and changes in the redox state in the photosynthetic electron transport chain are important indicators of operational retrograde signals. The function of operational retrograde signaling differs from biogenic signaling as it focuses on the adjustments and cellular homeostasis of the chloroplast in response to different perturbations and developmental cues [47].
Conversely, the most representative functional class identified in the chloroplast downregulated protein group (i.e., abundant in young leaves) were involved with cellular and organic biosynthetic processes, being the most representative in the "organic substance biosynthetic process" and the "nitrogen compound metabolic process". These results are under the developmental stage of the young leaves used here. As plastids are intimately related to developmental and environmental signals, chloroplasts from young leaves were expected to reflect the biogenesis processes, which involves large-scale gene transcription followed by synthesis, importing, and assemblies of biogenesis proteins, such as those involved in plastid translation apparatus and enzymes in the synthesis of isoprenoids and tetrapyrroles [48]. In turn, in mature leaves, redox-related proteins comprised most of the differentially regulated proteins; the only category suggested by gene ontology analysis in young leaves for "cell redox homeostasis" comprised only 3% of all downregulated proteins. A study on the developmental dynamics of maize leaves also found upregulated genes that encode enzymes involved in diverse biosynthetic processes in young leaves. Similarly, in mature leaves, the upregulated genes predominant encoded Calvin cycle, redox regulation, and the light reactions of photosynthesis [48].
For better insights into the changing patterns of the up and downregulated proteins, we classified the chloroplast proteins into three expression trends, as shown in Fig 3. The early responsive proteins comprised the subset with a drastic change in the abundance level in middle leaf tissues, remaining at similar levels in both mature and young leaves. Conversely, late responsive proteins showed a drastic abundance change only in the contrasting leaf tissue (young or mature leaves). We considered as almost linear responsive the proteins without drastic changes throughout the leaf development process, but with a significant change between mature and young tissues.
In conclusion, while investigating the chloroplast proteome changes over leaf development, we found significant alterations in the target proteome regarding the maturation process. Out of the 431 proteins isolated from the chloroplast material, 61 showed significant changes in the abundance throughout leaf development. There was no clear bias towards a unique change in the proteome composition; instead, we found distinct patterns of protein accumulation. In addition to pinpointing proteins with pronounced changes ranging the studied phenomenon, through a quantitative proteomics approach, we also specified the most particular biological and metabolic processes involving leaf development. However, as the proteome of any living entity is constantly changing, we also decided to investigate the undisturbed (i.e., without inducing anthropogenic alterations due to external stimulus) sub-cellular proteome dynamics by analyzing the seasonal changes occurring in the chloroplast proteome of Eucalyptus grandis plants cultivated at the State of Sao Paulo Forest Nursery between September 2019 and August 2020.

Transition seasons induce the most pronounced chloroplast proteome changes
Meteorological data revealed only small changes in the mean and maximum temperatures recorded throughout the study period, with the minimum diurnal temperature ranges observed during the summer season. Precipitation analysis revealed extremely dry fall and winter seasons and mean monthly rainfall of almost 10 mm over the summer (Fig 4A). Multivariate analysis of the E. grandis chloroplast proteomes suggested that both the fall and the spring-related samples presented the most dissimilar profiles of the studied seasons (Fig 4B). This is likely to be related to the phenological transitions commonly occurring in the species during such periods. In Populus species, phenology is largely related to seasonal changes in growth and development, such as bud flush in spring followed by vegetative growth in summer. In turn, during the fall, the trees presented growth cessation and bud formation, followed by leaf senescence and abscission, in addition to cold acclimation for the winter [32]. Our comparative proteome analysis revealed that 41 proteins were differentially regulated across the seasons of the year in E. grandis (S4 Table).
Hierarchical clustering analysis of DAP proteins suggests 5 different clusters according to their abundance patterns ( Fig 5A). As the transition seasons (fall and spring) presented the most dissimilar profiles, it is worth highlighting the clusters that showed the most expressive differences in these specific periods. Cluster 1 comprises 10 identifications with the lowest abundance values in the spring season. Proteins are related to light-harvesting complex, carbonyl detoxification, folding and modifications, carbohydrate biosynthesis, hydrolases, amino acid precursors, and ATP carriers. Conversely, Cluster 5 is composed of 6 proteins that showed higher abundance in the spring. Most of these identifications are associated with RuBisCO regulation to enhance photosynthetic capacity. Other identifications are responsible for lipid transport, protein modification, and translation ( Fig 5B).
Other studies have revealed that quantum efficiency seasonal pattern has minimal values in the spring, increasing over the summer, reaching maximal values in the fall, and finally decreasing in the winter [33,34]. Such scenario is demonstrated in Cluster 1, with proteins LHB1 and LHCB4, which participate in the first steps of photosynthesis and are very efficient in capturing light energy and transforming light into chemical energy [52]. Photosynthetic capacity in plants has also been reported to fluctuate with increased levels throughout leaf development, reaching its peak in the spring or early summer, and remaining stable or gradually declining over the summer [53][54][55]. Such pattern appears in Cluster 5, with RuBisCOrelated proteins (RCA, RBCS1, and RBCS2). Leaf fluorescence and photosynthesis are expected to increase over the spring, especially in chlorophyll content. In the fall, these parameters are reduced when associated with leaf physiological properties such as leaf aging and senescence [56]. Cluster 4 presented the highest abundance for the fall season (Fig 5B), with 7 identifications associated with the following processes: PSII regulation, chlorophyll biosynthesis, translation, phospholipid transport; in addition to two histones that play an important role in regulating transcription and DNA repair and replication.
As previously described, trees undergo growth cessation and establish dormancy before the onset of winter over a transition season. Low temperatures were found to facilitate photoinhibition but not necessarily to cause it [57]; in addition, downregulation of Calvin-Benson cycle activity and gas exchange precedes the reorganization of photosynthetic apparatus in thylakoid membrane during the fall [58]. All these events are partially related to PSII regulation and corroborate our analysis, where we found an upregulated protein related to the regulation of PSII in Cluster 4. Furthermore, a plant epigenome is very responsive to changes in the environment, and variations in epigenetic silencing mechanisms are among the key factors that assist plants to adapt to different climate conditions [59]. DNA methylation and histone modifications associated with chromatin remodeling were reported to be the major changes occurring throughout dormancy season transitions [60]. The up-regulation of two histones in the fall occurring in Cluster 4 could be associated with these epigenetics events.
Aiming to investigate which chloroplast proteins had the most pronounced synthesis in each season, we carried out a fold-change analysis, as shown in Fig 6. The red boxes indicate proteins with significant change equal or higher than 2-folds, blue boxes show significant change equal or higher than 5-folds, and white boxes point to proteins with change lower than 2-folds.
The most significant abundance increase had in both wet seasons (spring and summer) occurred in the RBCS1 protein. This is one of the small subunits of the RuBisCO enzyme, whose down-regulation has been observed in Cistus albidus [61] in drought periods, such as fall and winter. In Lolium perenne, the RbcS transcripts reached the highest levels over the spring, gradually decreasing between the summer and winter [62]. Similarly, all photosynthetic parameters measured in Quercus coccinea Muench (scarlet oak) showed a peak in the late summer, followed by a sharp decline in the fall [63].
In contrast, the most abundant proteins in both dry seasons (fall and winter) were MPH1, PORA, MEE32, TPI and AOR. MPH1 (proline-rich family protein) is a PSII-associated prolinerich protein that participates in the maintenance of normal PSII activity and protection of PSII from photooxidative damage under excessive light conditions [64]. The PORA protein (Protochlorophyllide oxidoreductase A) catalyzes the reduction of protochlorophyllide to chlorophyllide, subsequently producing chlorophyll a and b [65,66]. In addition to its role in chlorophyll biosynthesis, PORA can also exert a photoprotective action in the transition from dark to light [67], whereas the Triose phosphate isomerase enzyme (TPI) interconverts glyceraldehyde 3-phosphate (G3P) in the Calvin cycle placed in chloroplasts. Up-regulation of TPI in the winter has been reported in experiments performed on tubers and Eucalyptus grandis bark [68,69].
The bifunctional enzyme MEE32 (dehydroquinate dehydratase, putative/shikimate dehydrogenase) is responsible for steps three and four of the shikimate pathways, catalyzing the dehydration of dehydroquinate to dehydroshikimate and the reversible reduction of dehydroshikimate to shikimate, respectively. Transgenic lines of tobacco with suppressed MEE32 gene displayed severe growth retardation and reduced content of aromatic amino acids and downstream products such as lignin and chlorogenic acid [70]. In Eucalyptus camaldulensis, it confers aluminum tolerance [71]. AOR protein (Oxidoreductase, zinc-binding dehydrogenase family protein) detoxify α, β-unsaturated carbonyls by reducing a highly electrophilic α, βunsaturated bond using NAD(P)H [72,73]. It favors the detoxifying of stromal lipid peroxidederived reactive carbonyls (RCs) produced under oxidative stress in chloroplasts [73], in addition to protecting dark respiration and supporting plant growth during the night [74]. In conclusion, our findings suggest that transition seasons (spring and fall) induce the most pronounced chloroplast proteome changes over the year. Out of the 431 proteins isolated in the chloroplast proteome, 41 showed significant abundance changes in different seasons. We also found distinct patterns of protein abundance among the seasons and were able to pinpoint the most abundant proteins in both the dry and wet periods of the year.

The importance of studying chloroplasts metabolic changes
Chloroplast biogenesis starts with proplastids that are exposed to light and then transformed into mature chloroplasts, which differentiate from gerontoplasts during senescence. Leaf senescence can be triggered by both external and internal factors. Internal factors include age, reproductive growth, redox metabolism, and hormone levels [75][76][77]. Leaf premature senescence can cause productivity losses due to lower assimilatory capacity [78]. Investigating the complex molecular network behind the developmental processes is essential for understanding the continual transition of functional and regulatory metabolic pathways occurring multiple times throughout the perennials plants lifespan. It is also quite relevant for agriculture because [79,80] it may open new possibilities to increase crop yield and establish strategies to maximize biomass production.
An inter-organellar transcriptomic analysis of leaf senescence in Arabidopsis indicated that the chloroplast transcriptome exhibited the most pronounced changes throughout the plant lifespan [81]. These findings indicate that plastid compartment could be an important target to plant breeding programs that intend to delay leaf senescence and maximize biomass production. Regarding chloroplast metabolism, a large number of differentially expressed transcripts were identified in mature to senescence stages, and the biological processes associated with biogenesis at the growth to maturation stage showed higher expression, suggesting robust transcriptional coordination during early leaf development and more dynamic biochemical transitions during the senescence process [81]. At the protein level, we also observed higher abundance in biogenesis processes in young leaf tissue; in contrast, we detected a larger number of differential abundance proteins for the young stage rather than in mature tissues. This indicates that protein-level dynamics diverges from transcriptome dynamics. The plants also presented a strong translational coordination during early leaf development, and such biochemical networks seem more active during early stages.
Plants are sessile organisms and need to develop survival strategies due to changing environmental conditions. Perennial plants have been described to be strongly affected by seasonally varying availability of nutrients, water, light, and temperatures [30,31]. Seasonal climate variations can drive plant growth and development, favoring a metabolic adaptation of the life cycle with changes in the environment [29].
Our findings suggest that transition seasons (spring and fall) induce the most pronounced chloroplast proteome changes over the year in Eucalyptus grandis cultivated in tropical climate. For the fall, we found higher abundance of proteins related to epigenetic events that are very sensitive to external environment changes. Such epigenetic mechanisms are among the factors that help plants adapt to different environmental climates [59]. The combination of seasonal proteomics and epigenomics studies in natural environments can be important to reveal seasonal dynamics regulation at each level on plant responses plasticity.
Moreover, agrometeorological data revealed extremely dry autumn and winter seasons and mean monthly rainfall of almost 10 mm over the summer. Photosynthetic organisms have established several strategies to cope with dry environments, including short life cycle, enhanced water uptake and reduced water loss, osmotic adjustment, and antioxidant capacity [82]. Eucalyptus species have been reported to have developed a strategy involving lower photosynthetic activity during dry periods associated with resource reallocation through major changes in the gene expression of primary metabolism [83]. Identifying these metabolic changes allows to further investigate and understand the mechanisms behind plants adaptation and assist tree breeding programs.

Conclusions
Investigation of Eucalyptus grandis chloroplast proteome provided a comprehensive resource for understanding different processes underlying age-dependent and seasonal variation in an organellar level. Our chloroplast extraction and enrichment strategies were effective at reducing cross-contaminants abundance and allowed a deep investigation of the changes occurring in the chloroplast proteome throughout leaf development and the seasons of the year. Distinct patterns of protein accumulation were detected in both assays. The Eucalyptus grandis chloroplast proteome seems to increase metabolisms related to catabolic and redox processes as leaves age, while the accumulation of proteins involved in biosynthetic processes are downregulated. Most pronounced proteome changes occurred in the transition seasons of the year with a clear proteome bias towards photosynthesis over the wet seasons.
A deeper understanding on age-associated biological networks is important for developing strategies to maximize the biomass production. Furthermore, analyzing seasonal dynamics unravels the molecular mechanisms underlying metabolism plasticity and plant adaptation in spatiotemporal networks.   Table. Expression profile of the differentially regulated proteins identified in the Eucalyptus grandis chloroplasts during leaf development. Proteins were classified according to their abundance profile: "upregulated"-proteins with the highest abundance in mature leaves; "downregulated"-proteins with the lowest abundance in young leaves; "undefined"-proteins with undefined profile pattern.